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ABSTRACT 



In quantum mechanics textbooks the momentum operator is defined in the Cartesian coordinates and 
rarely the form of the momentum operator in spherical polar coordinates is discussed. Consequently 
one always generalizes the Cartesian prescription to other coordinates and falls in a trap. In this work 
we introduce the difficulties one faces when the question of the momentum operator in spherical polar 
coordinate comes. We have tried to point out most of the elementary quantum mechanical results, 
related to the momentum operator, which has coordinate dependence. We explicitly calculate the 
momentum expectation values in various bound states and show that the expectation value really 
turns out to be zero, a consequence of the fact that the momentum expectation value is real. We 
comment briefly on the status of the angular variables in quantum mechanics and the problems 
related in interpreting them as dynamical variables. At the end, we calculate the Heisenberg's 
equation of motion for the radial component of the momentum for the Hydrogen atom. 



Quantum mechanics is a treasure house of peculiar and interesting things. Elementary textbooks of quantum 
mechanics 0, [H, Q generally start with the postulates which are required to define the nature of the dynamical 
variables in the theory and their commutation relations. The choice of the dynamical variables is not that clear, as 
the coordinates in Cartesian system are all elevated to the status of operators where as time remains a parameter. 
More over in spherical polar coordinates only the radial component can be represented as an operator while the angles 
still remain as a problem. The difficulty of giving different status to the spatial coordinates and time is bypassed in 
quantum field theories where all the coordinates and time become parameters of the theory. But the problem with 
angles still remain a puzzle which requires to be understood in future. 

When we start to learn quantum mechanics, most of the time we begin with elementary calculations relating to the 
particle in a one dimensional infinite well, particle in a finite potential well, linear harmonic oscillator and so on. The 
main aim of these calculations is to solve the Schrodinger equation in the specific cases and find out the bound state 
energies and the energy eigenfunctions in coordinate space representation. While solving these problems we overlook 
the subtleties of other quantum mechanical objects as the definition of the momentum operator in various coordinates, 
the reality of its expectation value, etc.. In the last one or two decades there has been a number of studies regarding 
the self-adjointness of various operators [J. The aim of these studies has been to analyze the self-adjointness of 
various operators like momentum, Hamiltonian etc. and find out whether these operators are really self-adjoint in 
some interval of space where the theory is defined, if not then can there be any mathematical method by which we 
can make these operators to be self-adjoint in the specified intervals ? 

In the present work we deal with a much elementary concept in quantum mechanics related to the reality of the 
expectation values of the momentum operator, be it linear or angular. We do not analyze the self-adjointness of 
the operators which requires different mathematical techniques. To test the self-adjointness of an operator we have 
to see whether the operator is symmetric in a specific spatial interval and the functional domain of the operator 
and its adjoint are the same. In the article we always keep in touch with the recent findings from modern research 
about self-adjoint extensions but loosely we assume that the operators which we are dealing with are Hermitian. If 
something is in contrary we point it out in the main text. A considerable portion of our article deals with the analysis 
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of the fact that the expectation value of the momentum operator in various bound states are zero, a result which 
most of the textbooks only quote but never show. In the simpler cases the result can be shown by one or two lines 
of calculation, but in nontrivial potentials as the Morse potential, the Coulomb potential the result is established by 
using various properties of the special functions such as the associated Legendre and the associated Laguerre. 

The presentation of various materials in our article is done in the following way. Next section deals with the 
definition of the momentum operator and its properties. Section [Ml deals with the intricacies of the definition of the 
momentum operator in spherical polar coordinates and the problems we face when we try to mechanically implement 
the quantization condition, which is invariably always written in Cartesian coordinates in most of the textbooks on 
quantum mechanics. In section HVl we explicitly calculate the momentum expectation values in various potentials and 
show that in bound states we always get the expectation value of the linear momentum to be zero. Section IV1 gives 
a brief discussion on the Ehrenfest theorem when we are using it to find out the time derivative of the expectation 
value of the radial component of momentum in the case of the Hydrogen atom. We end with the concluding section 
which summarizes the findings in our article. 

Before going into the main discussion we would like to mention about the convention. We have deliberately put a 
hat over various symbols to show that they are operators in quantum mechanics. Some times this convention becomes 
tricky when we are dealing with angular variables as there the status of these variables is in question. The other 
symbols have their conventional meaning. As we are always using the coordinate representation sometimes we may 
omit the hat over the position operator as in this representation the position operator and its eigenvalues can be 
trivially interchanged. 



II. DEFINITION OF THE MOMENTUM OPERATOR AND THE REALITY OF ITS EXPECTATION 

VALUE 



From the Poisson bracket formalism of classical mechanics we can infer: 

[xi,pj] = ihSij , (2.1) 

where Sij = 1 when i = j and zero for all other cases, and i,j = 1,2,3. In the above equation Xi is the position 
operator and pj is the linear momentum operator in Cartesian coordinates. From the above equation we can also find 
the form of the momentum operator in position representation, which is: 

A = ■ (2-2) 

It is interesting to note that the above expression of the momentum operator also gives us the form of the generator 
of translations. This is because of the property: 

\p X) F(x)} = -ih^p-, (2.3) 

where F(x) is an arbitrary well defined function of x. The above equation ensures that the momentum operator 
generates translations along the x direction. 

Particularly in one-dimension the expression of the momentum operator becomes p x = —ih-^. We know that the 
expectation value of the momentum operator must be real. If we focus on one-dimensional systems to start with, 
where the system is specified by the wave-function ijj(x,t), the expectation value of any operator O is defined by: 

/oo 
tl>*(x,t)6tl>(x,t)dx, (2.4) 
-oo 

where ip* signifies complex conjugation of ip and the extent of the system is taken as — oo < x < oo. From the above 
equation we can write, 



(O)* = il>{x,t)0* il>{x,t)*dx. (2.5) 

J — OO 

If (O) = (O)* then the condition of the reality of the expectation value becomes: 

ip*(x,t)6ip{x,t)dx= / %j)(x,t)6* %j)(x,t)*dx. (2.6) 
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For a three-dimensional system the above condition becomes, 



t/>*(x, t)Oip(n, t)d 3 x 



t/>(x,i) O* V(x, t)*d 3 x. 



(2.7) 



Now if we take the specific case of the momentum operator in one-dimension we can explicitly show that its expectation 
value is real if the extent of the system is infinite and the wave- function vanishes at infinity. The proof is as follows. 
The expectation value of the momentum operator is: 



ip*(x, t) p x ip(x, t) dx — —ih I ip*(x,t) — ^— ^ — dx 



= —ih 



i)*(x,t)ip(x,t)\ 



dx 

oo 
— oo 



dii>(x,ty ' 

ip{x,t) dx 

OX 



(2.8) 



If the wave functions vanish at infinity then the first term on the second line on the right-hand side of the above 
equation drops and we have, 



ip*(x,t)p x ijj(x,t)dx 



-ih 



ih 



ip(x,t)— ^-T^r— dx , 
ip(x, t) p* tp(x, t)* dx . 



(2.9) 



A similar proof holds for the three-dimensional case where it is assumed that the wave-function vanishes at the 
boundary surface at infinity. 



III. THE EXPECTATION VALUE OF THE MOMENTUM OPERATOR IN CARTESIAN AND 

SPHERICAL POLAR COORDINATES 

In non-relativistic version of quantum mechanics we know that if we have a particle of mass m which is present in 
a time-independent potential we can separate the Schrodinger equation: 



ih 



at 



2 



h 



^(x,t), (3.1) 



into two equations, one is the time-dependent one which gives the trivial solution e n where E is the total energy 
of the particle, and the other equation is the time-independent Schrodinger equation: 

V 2 u(x) + ^-(E-V{x))u(x) = 0, (3.2) 

where u(x) is the solution of the time-independent Schrodinger equation and the complete solution of the Eq. (|3.ip 
is: 

ip(x,t) = it(x)e"^ ± . (3.3) 

In the case of the free-particle, where V(x) = 0, we have u(x, t) — e ik x where E — and k = |k|. The free- 
particle solution is an eigenfunction of the momentum operator with eigen value hk. Although if we try to find 
out the expectation value of the momentum operator as is done in the last section we will be in trouble as these 
wave-functions do not vanish at infinity, a typical property of free-particle solutions. But this problem is not related 
to the Hermiticity property of the momentum operator, it is related with the de-localized nature of the free-particle 
solution. 

In physics many times we require to solve a problem using curvilinear coordinate systems. The choice of our 
coordinate system depends upon the specific symmetry which we have at hand. Suppose we are working in spherical 
polar coordinates and the solution of Eq. (|3.2p can be separated into well behaved functions of r, 9 and 4> as, 



m(x) = u(r, 0, <f>) = R{r) Q{6) $(</>) 



(3.4) 
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If we try to follow the proof of the Hermiticity of the linear momenta components, as done in the last section, in 
spherical polar coordinates, then we should write: 



(p) = -ih J u*{r,9,4>)Vu(r,Q,4>)dT 
= -ih j i?*(r)9*(6»)$*(0) 



d eg d 



Or 



09 



r sin f 



R{r)0{9)<f>{4>)r 2 drdtt . 



(3.5) 



where in the above equation e r , eg, respectively are the unit vectors along r, 9 and <f>, and dCl — sin.6 d6 d<p. t is 
the volume over which we integrate the expression in the above equation. From the last equation we can write: 



(Pr) 



-ih 



|0(0)| 2 |<i>(0)| 2 dfi 



^ R * {r) im dr . 

dr 



(3.6) 



As 0(0) and <&(</>) are normalized, the integration: J n |0(#)| 2 |<I>(0)| 2 G?f2 = 1 and we can proceed as in Eq. (|2.9p as: 



("Pr, 



dR{r) 



-ih I r R*(r)- 
lo dr 



dr . 



r 2 R'ir)FHr)\;; - / ( 2rR* (r) + r 2 ^^ 



R{r) dr 



(3.7) 



If R(r) vanishes at infinity then the above equation reduces to, 



iPr) 



ih 



(Pr 



2n / 1 dR*(r) 

r z R(r) -r^- dr 



2ih 



dr 



r\R{r)\ 2 dr. 



2ih I r\R(r)\ 2 dr, 
'o 



(3.8) 



The above equation implies that (p r ) is not real in spherical polar coordinates. The solution of the above problem 
lies in redefining p r as is evident from Eq. (|3.8[) . and it was given by Dirac The redefined linear momentum 

operator along r can be: 



. d 1\ .1 d 
p r = -ih — + - = -ih- — r . 

or r I r or 



(3.9) 



This definition of the p r is suitable because in this form it satisfies the commutation relation as given in Eq. (|2.ip 
where now the operator conjugate to f is p r . The form of p r in Eq. Q3.9P shows that for any arbitrary function of r as 
F(r) we must still have Eq. (|2.3p satisfied. This implies that the modified form of p r is still a generator of translations 
along the r direction. Up to this point we were following what was said by Dirac regarding the status of the radial 
momentum operator. Still everything is not that smooth with the redefined operator as we can see that it turns out 
to be singular around r = 0, more over, although the radial momentum acts like a translation generator along r but 
near r — it cannot generate a translation towards the left as the interval ends there. 

In this regard we can state that the issue of the reality of the radial component of the momentum in spherical polar 
coordinates is a topic of modern research in theoretical physics 0, HI- It has been shown that the operator —ih-^ is 
not Hermitian and more over it can be shown [4| that such an operator cannot be self-adjoint in the interval [0, oo]. 
In some recent work [8] the author claims that there can be an unitary operator which connects —ih-S- to —ih--§- r, 

. . ...... or r or > 

and as the former operator does not have a self-adjoint extension in the semi-infinite interval so the latter is also not 
self-adjoint in the same interval. 

If we further try to find out whether (pg) and (p^) are real, then we will face difficulties. Working out naively if 
we claim that p^ = rs l n g -§g as suggested by the <f> component of Eq. (|3 . 5[) we will notice that 4>P<t> does not have the 



dimension of action. This means p^ or pg is not conjugate to <p or 9. This is a direct representation of the special 
coordinate dependence of the quantization condition. Only in Cartesian coordinates the variables conjugate to x, y 
and z are p x , p y and p z . Taking the clue from classical mechanics we know the proper dynamical variables conjugate 
to <j) and 9 are the angular momentum operators, namely Lg and L^. In general is given by: 



-ih 



d 



(3.10) 
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which can be shown to posses real expectation values by following a similar proof as is done in Eq. (|2.8p and Eq. (|2.9p . 
if we assume $(0) = <£>(27r). In this form it is tempting to say that we can have a relation of the form, 



[</>, = ih , 



(3.11) 



which looks algebraically correct. But the difficulty in writing such an equation is in the interpretation of <fi which 
has been elevated from an angular variable to a dynamical operator. In spherical polar coordinates both 9 and 4> 
are compact variables and consequently have their own subtleties. Much work is being done in trying to understand 
the status of angular variables and phases [13] , in this work we only present one example showing the difficulty of 
accepting <p as an operator. 

From the solution of the time-independent Schrodinger equation for an isotropic potential we will always have: 



$(0) 



JM<t> 



/2tt 



(3.12) 



where M = 0, ±1, ±2, •, •. Now if (f> is an operator we can find its expectation value, and it turns out to be: 



and the expectation value of (f> 2 is: 



Z7T 



2tt 



2tt 



2tt 



JAI4>„-iM<j> 



d<t>, 



b 2 e iM4 e -iM4 t 



(3.13) 



(3.14) 



Consequently Acp - 



2 ) — (cf>) 2 = Similarly calculating we get: 



to = t 

Z7T 

= Mh. 



2tt 



-iMcf) iM<fi 



d(j>, 



(3.15) 



as expected, and (L^) = M h . This implies AL^ = y (L 2 /) — (L^) 2 = 0. So we can immediately see that the 

Heisenberg uncertainty relation between 4> and L^, A^AL^ >h/2 breaks down. This fact makes life difficult and we 
have no means to eradicate this problem. 

Taking the clue from the <fi part we can propose that Lg is also of the form —ih-^g. With this definition of Lg let 
us try to prove its Hermitian nature as done in Eq. (|3.7p . Taking R(r) and $(0) in Eq. (|3.4[) separately normalized, 
we can write: 



= —ih 

= —ih 

= ih 

= (L e 



Q*(9)^± smOdd 
dO 

sm9Q*(d)0{6)\o- 



cos 9 0* (9) +sin6» 



dO*(9) 
d9 



ih 



d9 

cos9\@(9)\ 2 d9. 



ih / cos 9 1 6(0) | 2 d9 : 



Q{9) d9 



(3.16) 



The above equation shows that (Lg) is not real. The rest is similar to the analysis following Eq. (|3.8|) where now we 
have to redefine the angular momentum operator conjugate to 9 as [llj : 



Ln 



-ih 



d_ 

09 



■ cot ( 



(3.17) 
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Unlike the 4> case, Q(9) are not eigenfunctions of Lg. But the difficulties of establishing 9 as an operator still persists 
and in general 9 is not taken to be a dynamical operator in quantum mechanics. 

It is known that both 9 and <p are compact variables, i.e. they have a finite extent. But there is a difference between 
them. In spherical polar coordinates the range of <f> and 9 are not the same, < <f> < 2ir and < 9 < tt. This difference 
can have physical effects. As <f> runs over the whole angular range so the wave-function corresponding to it $(0) is 
periodic in nature while due to the range of 9, Q(9) need not be periodic. Consequently there can be a net angular 
momentum along the <p direction while there cannot be any net angular momentum along 9 direction. And this can 
be easily shown to be true. As the time-independent Schrodinger equation for an isotropic potential yields $(</>) as 
given in Eq. (|3 . 1 2|) similarly it is known that in such a potential the form of Q(9) is given by: 



9(6») = NgPfo (cos 0) 



(3.18) 



where Ng is a normalization constant depending on L, M and P^ I (cos9) is the associated Legendre function, which 
is real. In the above equation L and M are integers where L = 0, 1,2,3, •,• and M = 0, ±1, ±2, ±3, •, •. The 
quantum number M appearing in Eq. (|3.12[) and in Eq. (|3.18D are the same. This becomes evident when we solve 
the time-independent Schrodinger equation in spherical polar coordinates by the method of separation of variables. 
A requirement of the solution is — L < M < L. Now we can calculate the expectation value of Lg using the above 
wave- function and it is: 

(Lg) = —ifiNg J" Pt I (cOs9) ^M^Osg) + I cot QpL (cQS gfj gm m 

r dP L (co<*9) 1 f* 

= -ihNg 2 / Pm(cos9) M \ n J -sm0d0 + - Pm (cos 9)P^ (cos 9) cos 9 d9 

YJo df) 2 J 

(3.19) 



To evaluate the integrals on the right hand side of the above equation we can take x — cos 9 and then the expectation 
value becomes: 



(Lg) = -iTlNg 2 

1 



dP™(x) 



(1 - x 2 )* dx 



r M(x)^M(x) 



dx 

X 



-.dx 



(3.20) 



The second term in the right hand side of the above equation vanishes as the integrand is an odd function in the 
integration range. For the first integral we use the following recurrence relation [12l |: 



2 „<W 



(x<-iy- 



dx 



MxPt I (x)-(L + M)Pt I _ 1 (x). 



(3.21) 



the last integral can be written as, 

(Lg) = iUNg 2 



M I x(l-x 2 )~^P^(x)P^(x)dx 



As, 



(L + M) J (l-xt^M^M-iW^ 



P^(x) = (-l) L+M P^(-x), 



(3.22) 



(3.23) 



we can see immediately that both the integrands in the right hand side of the above equation is odd and consequently 
(Lg) = as expected. A similar analysis gives (L$) — MTi. It must be noted that the form of Lg still permits it to 
be the generator of rotations along the 9 direction. 

As the motion along <f> is closed so there can be a net flow of angular momentum along that direction but because 
the motion along 9 is not so, a net momentum along 9 direction will not conserve probability and consequently for 
probability conservation we must have expectation value of angular momentum along such a direction to be zero. 
In elementary quantum mechanics text books it is often loosely written that the solution of the time-independent 
Schrodinger equation is real when we are solving it for a real potential. But this statement is not correct. The reality 
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of the solution also depends upon the coordinate system used. Specially for compact periodic coordinates we can 
always have complex functions as solutions without breaking any laws of physics. 

Before leaving the discussion on angular variables in spherical polar coordinates we want to point out one simple 
thing which is interesting. In Cartesian coordinates when we deal with angular momentum we know that: 

[Li,Lj] = ie ijk L k , (3.24) 

where Li stands for L x , L y or L z . For this reason there cannot be any state which can be labelled by the quantum 
numbers of any two of the above angular momenta. But from the expressions of L^ and Lg we see that, 

[L <t> ,L e ]=Q, (3.25) 

and consequently in spherical polar coordinates we can have wave-function solutions of the Schrodinger equation 
which are simultaneous eigenfunctions of both and Lg as P^j (9) . 

For real V^x), we expect the solution of the time-independent Schrodinger equation it(x) to be real, when we 
are solving the problem in Cartesian coordinates. In all these cases the expectation value of the linear momentum 
operators must vanish. The reason is simple and can be understood in one-dimensional cases where with real u{x) we 
directly see that the integral u* (x) dU Q^ dx is real and so u* (x) p x u(x) dx becomes imaginary as p x contains 
i, as is evident from the first line in Eq. (|2.9|) . So if the expectation value of the momentum operator has to be 
real then the only outcome can be that for all those cases where we have a time-independent solution in a bounded 
region of space, with a real potential and working in Cartesian coordinates, the expectation value of the momentum 
operator must vanish. The above statement is true in curvilinear coordinates also, but in those cases the definition of 
the momentum operators have to be modified. This fact becomes clear when we write the relationship between the 
probability flux and the expectation value of the momentum operator. The probability flux for a particle of mass m 
is: 

j(x,t) = [f(x,t)V^x,i)-(Vf(x,t))^(x 1 t)] 1 



2m, 

Im(V>*(x,f)VV>(x,t)), (3.26) 



where 'Im' implies the imaginary part of some quantity. Most of the elementary quantum mechanics books then 
proceeds to show that: 



f dVj(x,t) = M (3.27) 
J m 



which is obtained from Eq. (|3.26|) by integrating both sides of it over the whole volume. From Eq. (|3.26[) we immediately 
see that if the solution of the time- independent Schrodinger equation is real we will have j (x, t) = and consequently 
from Eq. (|3.27[) . (p) = 0. But this statement is also coordinate dependent, which is rarely said in elementary textbooks 
of quantum mechanics. Eq. (|3.26|) evidently does not hold in spherical polar coordinates. If we take Eq. (|3.4|) as the 
solution in a general isotropic central potential and use the general form of V in spherical polar coordinates then it 
can be seen that j r (r, 0, </>, t) = for a real potential. But then Eq. (|3.27[) does not hold as here p r is simply the radial 
component of V and not as given in Eq. (|3.9p . and we know is not zero. The reason why Eq. (|3.26[) is not suitable 
in spherical polar coordinates is related to the fact that in deriving Eq. (|3 . 26[) one assumes that the probability density 
of finding the quantum state within position x and x + dx at time t is |^(x, t)\ 2 . But this statement is only true in 
Cartesian coordinates, in spherical polar coordinates the probability density of the system to be within a region r and 
r + dr, 8 and 8 + d6, 4> and (f> + d<f> is not \ip(r, 9, </>)| 2 but \ip(r, 9, </>)| 2 r 2 sin 5 and consequently the steps which follow 
leading to Eq. (I3.26|) in Cartesian coordinates are not valid in spherical polar coordinates. In general Eq. (|3.26|) will 
not be valid in any curvilinear coordinate system. 

The next section contains the actual calculations of the expectation values of the momentum operator in various 
cases where we have bound state solutions. In all the relevant cases discussed in this article it is seen that although 
(fix) = but (p x ) is not zero as it is related to the Hamiltonian operator. In all the cases we must have, 

((p x ) s ) = , s = odd integer . (3.28) 



The above equation can be guessed from the reality of the expectation value of the momentum operator. 
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IV. MOMENTUM EXPECTATION VALUES IN VARIOUS BOUND STATES 

In this section we will calculate the momentum expectation values in various bound states with stiff or slowly 
varying potentials. 

A. Particle in one-dimensional stiff potential wells 

1. Infinite square well potential 

In this case we consider a particle to be confined in region — h to 4 along the x-axis where the potential is specified 
by 

V(x) = oo, \x\ > —, 

= 0, \x\ < |. (4.1) 
In this case the solution of the time- independent Schrodinger equation, Eq. (|3.2|) , satisfies the boundary condition, 

. (-#)—(#)-». <"> 

and as the potential has parity symmetry about x = we have two sets of solutions, the odd solutions: 



2 . (2nirx\ 



u^(x) = Jy S m(^P), (4.3) 



and the even solutions: 

u^x) = ^ C os(^^). (4.4) 

In the above equations n is a positive integer. Both of these functions, u^\x) for the odd case and u„ \x) for the 
even case, are real and are not momentum eigenstates. But the momentum expectation values can be found out from 
the above solutions. For the odd solutions we have: 

/- \ ■* P (o)f Mn\x) J 

[Px) = -in / L u) n >{x) — — — dx, 



Aimrh f 2 . / 2nnx\ ( 2nicx 



L 2 J_l \ L J V L 

2 \ / \ 

= 0, (4.5) 

as expected. Similarly for the even solutions it is also easy to show that the expectation value of the momentum 
operator vanishes. 

2. Finite square well potential 

In this case, 

V(x) = 0, \x\>a, 

= -V ,\x\<a, (V o >0). (4.6) 

If we are not interested in the normalization constant of the bound state solution then the solution of the time- 
independent Schrodinger equation in this case is: 

u{x) ~ e~ K ^ , \x\ > a, 

r-j cos(kx) , \x\ < a, (even parity) 

~ sin(fcx) , |x|<a, (odd parity) , (4.7) 



9 



where, 



2m(-\E\+V ) 



2m\E\ 



(4.8) 
(4.9) 



In this case the expectation value of the momentum operator is: 

f + °° du(x) 
J-oo dx 

/ f—a poo 

e~ 2KX dx 



—ih 
0, 



^ x dx 



sin(fca;) cos(fca;) dx 



(4.10) 



where the first two lines of the above equation holds up to a constant arising from the normalization of the wave- 
function. In deriving the last equation we have taken the odd parity solution, but the result remains unaffected if we 
take the even parity solution as well. 



In this case the potential is: 



3. Dirac-delta potential 



V(x) = -V S(x), (V >0). 



(4.11) 



In this case there can be one bound state solution which is obtained after solving the Eq. (|3.2j) . Demanding that the 
solution u(x) satisfies the boundary conditions: 



u(x — — e) = u(x = +e) , 



du 
dx 



x—-\-e 



du 
dx 



-^^-u(x = 0) , 



where e is an infinitesimal quantity tending to zero, we get the form of the solution which is: 

u(x) = v^e , x < 0, 
= \fne KX , x > 0, 



(4.12) 
(4.13) 



(4.14) 
(4.15) 



where k 



_ mVg 



and the energy of the bound state is E = 



2/i- 



du(x) 

U(X) ; dX . 

dx 



The expectation value of the momentum operator in this case is: 

(fix) = -i^ 

= —ihn 
= 0. 



z 2kx dx 



-2kx 



dx 



(4.16) 



In this case, also from Hermiticity of the momentum operator we see that Eq. (|3.28|) holds true. 



B. Particle in one-dimensional slowly varying potentials 

1. Linear harmonic oscillator potential 

In the case of the linear harmonic oscillator we have: 

V(x) = -mw 2 i 2 , 



(4.17) 
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where u> is the angular frequency of the oscillator. The solution of Eq. (13. 2|) in this case, using the series solution 
method, yields: 

u n (q) = N n e-^H n (q), (4.18) 

where n = 0, 1,2, ■,■ and q — y/ax where a — H n (q) are Hermite polynomials of order n and N n is the 

normalization constant given by, 

^ = (^=^) 2 ■ (4-19) 



, ,/7rn!2 

The momentum expectation value in this case turns out to be, 
(px) = -iTi\fa [ u n (q) dU ^ q dq, 



= —ihyfaNj 



dq 

e-i 2 H n (q)^±d q - ^ q e-« 2 Hl{q)dq 
J-oo dq J_ 00 



= 0. (4.20) 

The first integral on the right side of the second line of the last equation vanishes because, dH ^ q ^ = 2niJ n _i(q) and 
consequently the integral transforms into the orthogonality condition of the Hermite polynomials. The second integral 
on the second line of the right side of the above equation vanishes because the integrand is an odd function of q. 

The linear harmonic oscillator (LHO) has some very interesting properties. To unravel them we have to digress a 
bit from the wave-mechanics approach which we have been following and follow the Dirac notation of bra and kets. 
The Hamiltonian of the LHO in one-dimension is: 

n 2 1 

H = + -muo 2 x 2 , (4.21) 
2m 2 

which can also be written as: 

H = ?i,bj I oJ a 4- 

2 

where a and d^ are the annihilation and the creation operators given by 



H = Kw (a)a+- ) , (4.22) 



if + - > a Wlr • ^ 

Zn \ moj J V in \ moj J 

It can be seen clearly from the above definitions that a is not an Hcrmitian operator. More over from the definition 
of the operators we see that, 

[a,fit] = l. ( 4 - 24 ) 

Conventionally the number operator is defined as: 

N^d^d, (4.25) 

and its eigen-basis are the number states |n) such that, 

N\n) =n\n). (4.26) 

The Hamiltonian of the LHO can be written in terms of the number operator and consequently the number states 
are energy eigenstates. In this basis the action of the annihilation and creation operators are as: 

a\n) = Vn|n- 1) , (4.27) 

d)\n) = Vn + 1 \n + 1) . (4.28) 

From the definitions of the annihilation and creation operators we can write the momentum operator as: 



Px 



i^_(o-at). (4.29) 
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From Eq. (|4.27[) . Eq. (|4.28[) and the above equation we can write the matrix elements of the momentum operator as: 

(n'\p x \n) = i^^^ (-Vn5 n >, n -i + \/n+ l8 n ',n+\) ■ (4.30) 

The above equation shows that the momentum operator can connect two different energy eigenstates. 

In the case of LHO, except the number operator states, we can have another state which is an eigenstate of the 
annihilation operator a. This state is conventionally called the coherent state and it is given as: 



9 00 

lot) = e 2 



where a is an arbitrary complex number. Now from Eq. (14.29|) we can find the momentum expectation value of the 
coherent state and it is, 



(Px) a = (a\p x \a) = Im(a) . (4.32) 

From the above equation we can see that although the expectation value of the momentum operator is zero in the 
energy eigen-basis but it is not so when we compute the momentum expectation value in the coherent state basis, 
which is essentially a superposition of energy eigenstates. It must be noted that the momentum expectation value is 
non zero only when the parameter a has an imaginary part. 

2. Poschl- Teller potential 

Among the potentials belonging to the hypergeometric class the Poschl- Teller potentials have been the most ex- 
tensively studied and used. This class of potentials consist of trigonometric as well as the hyperbolic type. The 
trigonometric versions have found applications in molecular and solid state physics and the hyperbolic variants have 
been used in various studies related to black hole perturbations. 

In the present work we use the trigonometric, symmetric Poschl- Teller potential given by: 

V[x) = Vbtan 2 (ax), (4.33) 

where Vq can be parameterized as: 

^o = ^ r A(A-l), (4.34) 

with for some positive number A > 1 and a is some scaling factor. The energy eigenvalues of the bound state solutions 
are: 

h 2 a 2 

E n = -- — (n 2 + 2nA + A), (4.35) 
2m 

and the solution of the time-independent Schrodinger equation is, 

u n (x) = N n ^cos(ax) (sin(ax)) , (4.36) 

where, 

1/2 



a(n + \)T(n + 2X) 

r(n + 1) 



(4.37) 



is the normalization constant and P£(x) is the associated Legendre function. At this point it is fair to point out that 
Ptf(x) is not the Legendre polynomial P^(x) appearing in Eq. (|3.18[) . as fi and v need not be integers as L and M. 
Pjf{x) is not a polynomial but the function appearing in the right hand side of Eq. (|4.36[) is a polynomial. 

Now as claimed in the text let us show that the momentum expectation value is indeed zero. Before we proceed 
let us simplify the notation a bit by calling /i — 1/2 — A and v = n + A — 1/2. Substituting z — ax we can write the 
momentum expectation value as: 

/it/2 j 
dz v/cos(z) P£ (sin(z)) — ( Vcos(^) P» (sin(z))) . (4.38) 
-tt/2 az V ' 
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Note the limits of the integration range from ir/2 to — ir/2 since at this value the potential becomes infinity hence we 
need not consider the integration range to be the whole real line. For the sake of convenience let us make a change 
of variable; letting y = sin(z) the above integral becomes: 



(p x ) = -ihNl 



+1 dy (1 - tff/ipj;(y)± [(1 - y^P^y) 
-i ay L 



Taking the derivative inside the integral we get: 

H-l 



(Pa 



iUNt 



dy 



It is known that for associated Legendre functions \lc 



(4.39) 



(4.40) 



P£(-x) = cos[{p + i/)tt] P£(s) - - wo[(jjl + v)ir] Q$(x) . 



(4.41) 



where Q^(x) is the other linearly independent solution of the associated Legendre differential equation. As in our 
case fi + v = n so P{f(x) will have definite parity. As Pj?(x) has definite parity so the contribution of the second term 
in the above integral vanishes since the total integrand is an odd function. The first integral is similar to the one in 
Eq. (|3.20p and, due to the typical parity property of Py{x) as shown in Eq. (|4.41ll . it also vanishes. Consequently we 
have {px) = as expected. 



3. Morse potential 



Diatomic molecule is an exactly solvable system, if one neglects the molecular rotation. The most convenient model 
to describe the system, is the Morse potential [bH ]: 



V(x) = D(e 



-20x 



- 2e 



(4.42) 



where x = r/ro — 1, which is the distance from the equilibrium position scaled by the equilibrium value of the 
inter- nuclear distance r . I? is the depth of the potential, called dissociation energy of the molecule and f3 being a 
parameter which controls the width of the potential. 

In terms of the above scaled variable x, the time-independent Schrodinger equation becomes: 



fi 2 d 2 u(x) 



D(e 



-2/3 x 



2e~ l3x )u(x) = Eu{x) 



(4.43) 



2fj.ro dx 2 

Here /i is the reduced mass of the molecule and the corresponding bound state eigen function comes out to be: 

u£(£) = Ne-^ s/2 L s n (0 , (4.44) 



where the variables are described as, 



and 



Z = 2\e~ v ; y = (3x; < £ < oo , 



= 0,1,..., [A -1/2], 



(4.45) 



(4.46) 



which is nothing but the quantum number of the vibrational bound states. Here [p] denotes the largest integer smaller 
than p, thus total number of bound states is [A — 1/2] + 1. The parameters, 



' 2p,Drl 



I3 2 Ti 



and s = 



'(3 2 h 2 



E . 



(4.47) 



satisfy the constraint condition s + 2n = 2 A — 1. We note that the parameter A is potential dependent and s is related 
to energy E. In Eq. (|4.44[) , L s n (y) is the associated Laguerre polynomial and N is the normalization constant (l5| : 



N 



f3(2\ - 2n - l)T(n + 1) 
r(2A - n)r Q 



1/2 



(4.48) 
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We are looking for the expectation value of linear momentum for a vibrating diatomic molecule, and its expression 



is: 



(Px) = -ih I Un(£)—u n (£)dx. 

-OO ®% 



(4.49) 



In terms of the changed variable £ — 2Ae " x the integration limit changes to oo to and the expectation value 
becomes: 



(p x ) = -ih I K(Q—v, n ((;)d£ 



= ihN 



i r s r 

J\-ez'L> n (z)±L' n (£)dt 



ihN z 



Integral I\ is the orthogonality relation of the associated Laguerre polynomials, which is: 

T(s + n + 1) 



Jo 



r(n+l) 



(4.50) 



(4.51) 



To evaluate the second integral one uses the normalization integral of Morse eigenstates. The normalization relation 
is: 



\N\ 2 r 



e-^- 1 (L;(0) a de = l. 



(4.52) 



The above integral involving £, is explicitly l<i. N, being the normalization constant as given in Eq. 14.481 Thus it is 
very straight forward to evaluate J 2 from the above relation as, 



h = 



T(n + s + 1) 



sT(n + l) 

The last integrand 73 includes a differentiation which can be written as [l| 
Writing the right hand side of the above equation as a summation [T3] : 

n 
m— 

and substituting the derivative term in integral I3 we obtain: 

n—l 



m=0 ^° 



(4.53) 



(4.54) 



(4.55) 



(4.56) 



In the above integral m 7^ n because m can go only upto (n — 1). Thus the integral vanishes. Now let us see what is 
the expectation value of momentum observable, after evaluating the three integrals above. Substituting the non-zero 
values I\ and I2 in Eq. 14.501 it is clear that the expectation value of momentum is zero as has been expected. 



C. Position expectation values for various potentials 



After a thorough discussion about the momentum expectation values for various solvable one-dimensional potentials, 
it is worth spending some time discussing about the average position of the particle inside the bound states. Among 
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all the above examples, in each case we had V(x) — V(—x) except the Morse potential as Morse potential is not an 
example of a symmetric potential: V(x) ^ V(—x). 

In deriving the expectation values of momentum for above symmetric cases, we often considered that the integrals of 
odd functions over the symmetric limits vanishes. This result does not hold true for the asymmetric Morse potential. 
Already we have shown that the momentum expectation value: < p >= for all the above potentials. When it 
comes to the expectation values of position, one can easily see that < x >= for symmetric potentials whose centers 
are at the origin. On the other hand if this is not the case, suppose the infinite square well is defined in the range 
< x < L also then the expectation value of position does not vanish. It becomes L/2. Thus, more accurately the 
average position of the particle is dependent on the symmetry of the potential where as the average momentum is 
solely guided by the reality of it's eigenvalues and consequently it is zero always. 

Below we will briefly discuss how the asymmetry of the potential affects the expectation value of x in the case of 
the Morse potential. The expectation value of the position operator is: 



The eigen function and the variables are respectively substituted from Eq. (|4.44|) and Eq. (|4.45|) . We obtain 
... N 2 r 



ln(2A) 



(4.57) 



(4.58) 



The first integral is already been obtained in Eq. (|4.53|) . This result is independent of the quantum number n. The 
second integral (say /) is not that straight forward, because it contains associated Laguerre polynomial, logarithm, 
exponential and monomial functions. Here at best we can evaluate the integral atleast for some specific n as, n = 
or n = 1, when the Laguerre polynomial is respectively replaced by 1 and (— £ + s + 1). For the ground state wave 
function (n = 0), I would be 



I n =o — 



which can be written in terms of ^(s) and T function [If 



n=0 



(4.59) 



(4.60) 



where, ^(s) is the logarithmic factorial function, defined as d ( ln ( s ) ! 



- = \I/(s). For n — 0, first integral reduces 
to r(s) from Eq. (|4.53[) . Above two evaluations gives the ground state expectation value: 



\X)n=0 



1 



[ln(s + 1) - *(«)] 



For n = 1, one can proceed in the same way 

= ~P 

N 2 



(x} 7 



e-«f +1 ln(£K + (s + lj 



poo poo 

/ e-e? - 1 ln(£)d£2( S + 1) / e"«r 
Jo Jo 



In(£K 



(4.61) 



(4.62) 



[r(s + 2)*(a + 2) + (s + l) 2 r(s)*(s) - 2(s + l)T{s + l)*(s + 1)] , 



which simplifies to give the expectation value corresponding to the second eigen state: 



{X) n =l 



1 

^ol 



ln(s + 3) - *(s + 2) + 



(8 + 2) 



(4.63) 



Other expectation values for n > 1 can also be obtained in a similar fashion. 

The important point which is to be noted here is, though the average momentum vanishes, the average position is 
non-zero for Morse potential and remain so, irrespective of the choice of coordinate origin. This result is also true for 
all eigen states of the same Hamiltonian. 



D. Momentum expectation value for a three-dimensional slowly varying spherically symmetric potential 

In three dimensions, for a spherically symmetric potential the solution of the Schrodinger equation is given in 
Eq. (|3.4j) . Here we have assumed that the variables can be separated. The expectation values of Lg and have been 
evaluated in section IIII1 In this section we take the case of the Hydrogen atom and calculate the expectation value of 
the radial component of the linear momentum. 
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1. The Hydrogen atom 



In this case, 



V{f) 



where e is the electronic charge and r = \/ x 2 + y 2 + z 2 . Now we have to write Eq. 
and the solution of the time-independent Schrodinger equation is: 



inLM(r,0,<fi) = N r R nL (r)Y LM (0,< 



N r e ~ r/nao 



2r 



C 



2L + 1 
i-L-1 



( 2r 
\na 



Y, 



LM( 



(4.64) 

in spherical polar coordinates 
0) , (4.65) 



where eto = is the Bohr radius and m is the reduced mass of the system comprising of the proton and the electron. 
n is the principal quantum number which is a positive integer, C^_^_i{x) are the associated Laguerre polynomials, 
Yl m(9, 4>) are the spherical-harmonics, and N r is the normalization arising from the radial part of the eigenfunction. 
The values which L and M can take is discussed in section 11111 The radial normalization constant is given by: 



N r = 



( 2 \ 3 (n-L-l)\ 
\nao J (n + L)\2n 



1/2 



The spherical-harmonics are given by, 

Y LM (t 



2L + 1 (L-M)l 



-1 1/2 



4tt (L + M)\ 



Pfc(cos8)e 



iM<j> 



(4.66) 



(4.67) 



where Pff{cosO) are the associated Legendre functions. It is noted that although the Coulomb potential is a real 
potential but the solution in spherical polar coordinates is not real, e tMc ^, is complex. The spherical-harmonics are 
ortho-normalized according to the relation, 



d6d(t> sin 9Y LM {9, <£)F £jt (0,0) =& L L S . 



=0 J<t>=o 



J LL AIM ' 



(4.68) 



Let us write the eigenf unctions in terms of dimcnsionless quantity: p = 2r/nao = ar. Also we define k = (2L + 1) 
and n r = (n — L — 1) for the sake of convenience. With this amount of notational machinery the eigenfunctions can 
be written as: 



u n l m (r, 9, 4>) = N r R nL (p) Y LM (9, < 



(4.69) 



The radial momentum expectation value in this case is not given by —ih(J^), its form is (already discussed in section 



(p p ) = - l hN 2 dpp 2 R* nL {p) 



— + - 

dp p 



Rul(p) / dtl[Y LM {9A)? 



(4.70) 



Where TV 2 = N 2 /a 2 . The integral for the spherical harmonics yields identity. The radial expectation value then 
becomes, 



r°° ( 1 d 

(P P ) = - ihN2 J o dp! [ --e-Pp k +'l£ k nr (p)} 2 + (L + l)e-Pp k [£l(p)] 2 + e^p k+1 Cl(p)-[Cl(p)} 



Using the recurrence relation [If 



dp 



the expectation value integral acquires the form: 

(p p ) = -ihN 2 (°° dp Lle-Pp k+1 [C k n (p)] 2 + (n r + L + l) e~" p k [C k n (p)f + e^p k L k {p)C k {p) 



(4.71) 



(4.72) 



(4.73) 
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The third contribution of the becomes zero from the orthogonality property of the associated Laguerre polynomials 
as given in Eq. (|4.51[) . The contribution from the second term can also be found similarly. To find the share of the 
first term we make use of (l9j |: 

l°° dpe-"p k+ l [C k rir ( P )f = (2n r + k + 1). (4.74) 

Jo n r \ 

Collecting all the contributions we get the radial expectation value to be zero as expected. 



V. A DISCUSSION ON HEISENBERG'S EQUATION OF MOTION AND EHRENFEST THEOREM 

The time evolution of any operator O in the Heisenberg picture is given by: 

where H is the Hamiltonian of the system. The Hamiltonian of a quantum system comprising of a particle of mass 
m is given by: 

H=^ + V(x). (5.2) 
2m 

From the above two equations we can write the time evolution of the momentum operator in one dimension, in 
Cartesian coordinates as: 

% = ±\p x ,H] = -±V(x), (5.3) 
at in ax 

which is the operator version of Newton's second law in a time independent potential. Now if we take the expectation 
values of both sides of Eq. (|5.3[) in any basis we get: 

^~<5"<*>>- M 

and historically the above equation is called the Ehrenfest theorem, which was deduced in a different way by P. 
Ehrenfest. Using the Ehrenfest theorem we can deduce that the rate of change of the expectation value of the 
momentum operator is zero in the case of the linear harmonic oscillator. In the case of the linear harmonic oscillator 
we have: 

jLytx) — mu 2 x , (5.5) 
ax 

and it can be trivially shown that (x) = 0. This directly implies that, 

d(p x 



dt 



0, (5.6) 



for the linear harmonic oscillator. The above equation shows that the expectation value of the momentum along x 
direction is constant, and this constant is zero is known from other sources. 

Next we focus on the Hydrogen atom. The Hamiltonian of the Hydrogen atom is: 

S 2 19 2 1 e 2 

2m r or z 2mr z r 

where, 

\ sm oO oO sin 8 ocp z ) 

whose eigenvalues are of the form h 2 L(L + 1) in the basis Y^m(S, 4>)- In the expression of the Hamiltonian m is the 
reduced mass of the system comprising of the proton and electron. Next we try to apply Heisenberg's equation to the 



17 



radial momentum operator. Noting that the first term of the Hamiltonian is nothing but p 2 the Heisenberg equation 
is: 



dp r 



\2_ 

2m 
L 2 



ld_ 

r dr 



1 



1 d 1 

r dr r 



(5.9) 



The above equation is the operator form of Newton's second law in spherical polar coordinates. Next we evaluate the 
expectation value of both the sides of the above equation using the wave-functions given in Eq. (|4.65p . We know, 



1 



1 



3 2 



a o( L - 



7 3 7-, 3 



L(L + i)(L + l) 



(5.10) 
(5.11) 



Using the above expectation values in Eq. (|5.9p and noting that (L 2 ) = h 2 L(L + 1) we see that the time derivative of 
the expectation value of the radial momentum operator of the Hydrogen atom vanishes. 

The above analysis shows that the form of the Ehrenfest theorem as given in Eq. (|5.4p is only valid in Cartesian 
coordinates. In the case of the Hydrogen atom if we used Eq. (|5.4[) we should have never got the correct result. 



VI. CONCLUSION 



In the present work we have emphasized on the reality of the momentum expectation value and using the reality 
of the expectation value as a bench mark we did find out the form of the momentum operator in spherical polar 
coordinate system. We found that most of the concepts which define the momentum operator in Cartesian coordinates 
do not hold good in spherical polar coordinates and in general in any other coordinate system. The reason being 
that whenever we do an integration in curvilinear coordinates the Jacobian of the coordinate transformation matrix 
comes inside the picture and the Cartesian results start to falter if we do not change the rules appropriately. The 
forms of the momentum along the radial direction and the form of the angular momentum operators are derived in 
section Hill The status of the angular variables was briefly discussed in the same section. We explicitly calculated the 
expectation values of the momentum operator in various important cases and showed that the expectation value of 
the momentum operator do really come out to be zero as expected. Although the expectation value of the momentum 
operator vanishes in most of the bound states, with a real potential, the expectation value of the position is not 
required to vanish. The expectation value of the position operator is directly related with the parity property of the 
potential which was briefly discussed in subsection IIV CI At the end we calculated the Heisenberg equation of motion 
for the radial momentum operator for the Hydrogen atom and showed its formal semblance with Newton's second 
law. It was also shown that if we properly write the Heisenberg equation of motion in spherical polar coordinates 
then Ehrenfest 's theorem follows naturally. 

In short we conclude by saying: 

1. the forms of the various momentum operators, in most of the coordinate systems, in quantum mechanics can be 
obtained by imposing the condition of the reality of their eigenvalues. The form of the probability conservation 
equation and Ehrenfest theorem must be modified in curvilinear coordinates to yield meaningful results. 

2. There are obvious problems in elevating the status of angular variables to dynamical variables in quantum 
mechanics. 

3. For compact variables, if the variable is periodic the expectation value of the angular momentum conjugate to 
it is non-zero. If the compact variable is not periodic then the angular momentum conjugate to it must vanish. 

4. The momentum expectation values in cases of bound state motions vanish, whereas the position expectation 
values in those cases depends on the symmetry of the potential. 
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